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The Method of Repeated Application of a Quadrature Formula of 
Trapezoids and Rectangles to Determine the Values of Multiple Integrals 



Abstract: The work deals with the construction of multi-dimensional quadrature formulas based on the method of repeated 
application of the quadrature formulas of rectangles and trapezoids, to calculate the multiple integrals' values. We've 
proved the correctness of the quadrature formulas. 

Keywords: multi-dimensional integrals, etthaliire formulas, re-use method. 



The paper is devoted to developing the multiple cubature formulas for the values of n-fold integrals calculation by 
means of methods of repeated application of, quadrature formulas of trapezoid and rectangles. 

According to mathematical analysis multiple integrals can be calculated by recalculating single integrals. 
The essence of this approach is the following. Let the domain of integration is limited to single-valued continuous curves [1] 
y = q>(x), y = y/(x) (q>(x) < q>(x)) 
and two vertical lines X = a, X = b (Fig.l). 
Placing the known rules in the double integral 
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I. INTRODUCNION 
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limits of integration, we obtain 




(2) 



Then 




(3) 



y = q>(x 




y = (p(x 
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Fig.l 



Applying to a single integral in the right-hand side of (3), one of the quadrature formulas, we have 




where - X ; e [a,b] (i = 1,2,..., ri) and A ; are some constant coefficients. In turn, the value 
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F(x ( )= f f(x i ,y)dy 



can also be found on some quadrature formulas 



F(x i ) = f j B ij f(x i ,y J ), 



where 5- iy - appropriate constants. 
From (4) we obtain 



\\f(x,y)dxdy 



XXA.5,/(x„y.), 



(5) 



where A j and B {] iy - are known constants. 

Geometrically, this method is equivalent to the calculation of the volume I, given by the integral (2) by means of 
cross-sections. 

For cubature formulas of the type (3), the general comments, relating to the computation of simple integrals, are 
valid with appropriate modifications. 
In three dimensions (3) has the form 

JJJ /(*> y. z)dxdydz = ZZZ A^.C,^/^ ,y p z k ). 

{D) <=1 M k=l 

Now these results are generalized to compute the values of multiple integrals 



Where A-^ A.^ ^ - are known constants. 

This approach can be applied to calculate the approximate value of multiple integrals using appropriate quadrature 
formulas. In this case, we apply the quadrature formula of trapezoid. 



Is continuous in a bounded domain D. Then the formula for the approximate values of the n-fold integral calculating 




II. ALGORITHM DISCRIPTION 



Theorem 1 . Let the function 



y = f(x l ,x 2 ,...,x„) 




(6) based on the trapezoidal 



rule has the form 




Here D is n-dimensional domain of integration of the form 

a,. < x ( . <A i , i = l,n; (7) 




and it is believed that each interval (5), is respectively, divided into n l , n 2 ,..., tl n parts. 

Proof. Approximate formula for calculating the values of the one -dimensional integral in this case is 

/, = ]f(x)dx= h ^[f(x 0 ) + f(x l )]= h ±Yf i , (8) 




where, h x = , n - the number of the interval divisions into n parts. In (4) n = 1 

n 

Now we define the approximate calculation formula for double integrals 




(9) 
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where D is the two-dimensional area of integration. 

It is generally believed that the region of integration in this case is within the rectangle, that is, 

a<x<b, c<y<d. 

Then we rewrite the integral (7) in the form 

J 2 =\dx\f(x,y)dy. (10) 

Using the trapezoidal rule to calculate the inner integral, we obtain 

_ h y 
2 ~~2~[ 

Now to each integral we use the trapezoid formulas and the result is 

Ji = \ W ' y ° } + /(Xi ' y ° )] + h } + f(Xi ,yi 4 = 

' (12) 



j/(x, y)dx + \f(x, y)dx 



jfcjo)+/^.Jo)+MJi)+MJi)]=jEE/(^^' 



Similarly, for the approximate calculation of the values of the triple integral we obtain the following formula: 

+ f(x 0 , y, , Zl ) + f( Xl , y Q ,z 0 ) + fix, , y 0 , Zl ) + f{x x ,y 1 Z 0 ) + f(x 1 , y lZo ) + f(x, , y lf zj] = (13) 
h h h 1 1 1 

O , =0 7 =o k=Q 

C-c 
where h z = . 

The four-dimensional integral is defined with 16 members: 

J, = J J J J f{x,y,z, u)dxdydzdu = h * h & h * [f(x 0 ,y 0 ,z a ,u 0 ) + f(x 0 ,y a ,z a ,u l ) + f(x 0 ,y 0 ,z i ,u 0 ) + 

+ f(x 0 ,y 0 ,z } ,u 1 ) + f( XQ , yi ,z 0 ,u 0 ) + f(x 0 ,y l ,z 0 ,u ] ) + f(x 0 ,y l ,z 1 ,u 0 ) + f(x 0 ,y ] ,z l ,u l )+ (14) 

+ f(x l ,y 0 ,z 0 ,u 0 ) + f(x l ,y 0 ,z Q ,u 1 ) + f(x 1 ,y 0 ,z u u 0 ) + f(x l ,y 0 ,z 1 ,u 1 ) + f(x l ,y 1 ,z 0 ,u 0 ) + 

n hhhh I''' 
+ f(x„y 1 ,z 0 ,u l ) + f(x„y 1 ,z l ,u l )]= ^ ^ ZZZZ/fr/.y^*."/); 



Considering (8.10), we present formulas for the approximate calculation of the n-fold integral on the trapezoidal: 

/„= Jl|..j/(x 1 ,x 2 ,...,xj&A 2 ^3"A = ^n^zzi"i/K' x ^'"-' jc »i) (i5) 

Then for the formulas (6) and (8) - (11) we derive the general trapezoid formula. For this purpose, one-dimensional integral 

for integral intervals [a, b] is divided by n equal parts: [x 0 ,X 1 ],...,[x n _ 1 ,X n ], X i =X Q +ih x . 

The two-dimensional integral of the [a, b] x [c, d] region of integration is divided, respectively, by n and m parts: 

[x 0 , x r ] ,. .. , [x n _ x ,xj; [y 0 ,yj,...,[ y m _ y ,y m ] 

x = x 0 + ih , 

(16) 

yj = y 0 +J h y- 

In the same way ,the three-, four-, etc. n-dimensional integration region is divided by the corresponding parts. 
In this approach, the formula (6) takes the form 
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(17) 



. /o+/„ 



The general trapezoid formula for double integrals calculating has the form 

J 2 = £ £ k; + + f i+ u + Wi 1 = 

4 i=0 ;=o (1 8) 

The following is a general formula for the trapezoidal three-and four-fold integrals calculating 

J s =]]]fix.y,z)dxdydz = ^^figgk,,* < /.,.<•: ' Aj^ + 

abc d h h h h v 

J 4 =\\\\f(x,y,z,u)dxdydzdu =^^[f u , k , l + f iJXM + f iJMU + 

+ fiJMl,M + fi.j+Ul + fi,j +W+ l + AniMU + fiJ + lMl,M + f M J,k,l + f M ,m + 

1 (20) 

+ f M ,jMU + f M JMl,M + f M ,j +W + //+ w+ u./+i + Wu + u ++// + u + u + i.mj = 
h h h h "- 1 m ~ x p- 1 (-iiiii 

J-U i=0 y = o a = o /=(),■=()/, =o/ 3 =0 i 4 =0 

Based on the above formula, the approximate calculation of n-fold integral, using the formula (D), takes the form 

z /=1 p 1 =Op 2 =Qp 3 =0 p„ =()/,-()/,-() /„-() 

Thus the theorem is completely proved. 

Theorem 2. Approximation formulas for determining the values of the double integral, obtained by the repeated use 
of the trapezoid quadrature formula is calculated as follows: 



JJ fix, y)dxdy = ^pZSX/o 



where 



1. 2 4 ... 2 1 
2 4 4 ... 4 2 
2 4 4 ... 4 2 



2 4 4 ... 4 2 
12 2 ... 2 1 

Proof. Let's use the formula (12). Revealing this amount and similar terms, we obtain the following formula: 



(22) 



i23i 



h.li, 



XXf/., •./;,, '/.,... •/..,,,! 



= [/> 1 + 2/ 12 + 2/ 13 + ... + 2f ln _, + f l n 

2/21 + 4 /22 + 4 /23 + ••• + 4/ 2 „_! + 4f 2 „ 



4/ rtl +4/„_ 1>2 +4/„_ 1>3 +... + 4/„_ 1 „_ I +4/„_ 1 „. 

/«.! + 2/„,2 + /„.3 + - + 2/„,„_ 1 + /„.„ ] 

The (17) can form the matrix elements - A.^ in the form (16) and the result (15), indicating the proof of the theorem. 
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Theorem 3. Approximate formula for calculating the value of triple integrals, obtained by the repeated use of the 
trapezoid quadrature formula, is calculated as follows: 

hhh » » s 
x,y,z)dxdydz = y ZZZWp - (25) 



where 



12 2 2 
2 4 4 4 
2 4 4 4 



2 4 4 4 

12 2 2 

2 4 4 ... 

4 8 8 ... 



4 8 8 
2 4 4 



... 2 1 
..4 2 
..4 2 



..4 2 
... 2 1 



4 4 2 



Proof. Let's use (13). Revealing this amount and similar terms, we obtain 
./', , + 2/,, , + 2f 1X3 + 2f l l 4 + ... + 2J\ , „ , + y, „ 
2/ w + 4 /.,2.2 +4/ 1A3 +4/i A 4+... + 4V,+2/ 1A , 



2/i,„-i.i +4/,,„_ 1 , 2 +4/ ln _ 13 +4/ 1>n _ 1>4 + ... + 4/, „_,„_, + 2/ 1>n _ 1 „ 
/w + 2/i,„,2 + 2/^.3 + 2/, „ 4 + ... + 2f X n n _ x + f l n n 
2 / 2.1.1 +4/ 2 .,, 2 +4/ 213 +4/ 214 + ... + 4/ 21 „_ 1 + 2f 2Xn 

4/2,2,1 +8/2.2,2 +8/2,2,3 +8/2,2.4 +- + 8/ 2 .2.„-l +4/2,2,,, 



4 /2.,o-,i +8/2, 10 -,2 +8/2,10.3 +8/ 2 , 10 -, 4 +- + 8 / 2 ,io-, B -i +4 / 2 , 

2/2,11,1 +4/2.11.2 +4/2.11,3 + 4 /2,l V* +- + 4/2,1 + 2/2.1 1.„ 



2 /n-i,i,i +4/„_ w>2 +4/„_i.,,3 +4/„_ 1>M +- + 4/„_ WiB _ 1 + 2f n _ 1Xn 
4/„-i,2,i +8/„_,,2, 2 +8/„_i, 2 ,3 +8/„_ 124 + ... + 8/„_ 12 „_ 1 +4/„_ 12 „ 



4/„_ 1 ,„_ 1 ,i + 8/ n _ lin _ li2 + 8/„_, .„_, , 3 +8/„_,„_ 14 + ... + 8/,,,! .„_,,„_, +4/„_ 1 „_ 1 „ 
2/„„.„, 1 +4/„_ 1 ,„, 2 +4/„_,, n , n +4/ n _ 1 ^ + ... + 4/„_ 1 „„_ 1 +2/ n _ lBn 

/,,u + 2/„, 1>2 + 2/ n>1>3 + 2f nX4 + ... + 2/„ ,„_, + f nX „ (27) 
2/„.2.1 +4/„. 2 . 2 +4/„. 2 .3 +4/ n>2 , 4 +- + 4/„ j2 ,„-i +2/„. 2 .„ 



2/„,„_ u + 4/„.„_ 1 , 2 + 4/„ , n _ 1>3 + 4/„_ n _ lj4 + ... + 4f n + ... + 2f n ^ n 

fn,n,! + 2/ B ,„,2 + 2/„, K , 3 + 2/„ ,„ >4 + ... + 4/„ „ + 

From (19) we form the A ijk elements, resulting in the relation (18), which shows the proof of the theorem. 
Theorem 4. Approximation formulas for the values determining of the quadruple integrals, obtained by the repeated use of 
the trapezoid quadrature formula, is calculated as follows: 

rrr r h r hhh "4"^ 

\\\\f{x,y,z,u)dxdydzdu = * * ZZZIX^/ , (28) 

( D ) it) , =u j =Q k=0 ; =0 
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A n _ 


.1 A.-1.2 ■ 


■ A,-i,„- 


i A-i.„ 


A,., 


A,2 


■ A.„-i 


An 



1 = An = An = 



A, = An = 



2 4 



8 16 
8 16 



Proof. Disclose the amount (14) and obtain the following relations: 
/iin + 2/1112+ 2/1113 + 2/iii4 + - + 2/iii^-i+/iiv. 
2/nzi + 4/, 122 + 4/, [23 + 4f x 124 + ... + 4/, u „_, + 2/, u „ 
2/1131 + 4/ 1132 + 4/ 1I33 + 4/ I134 + ...+ 4/j Un _j + 2/ m „ 



16 8 
16 8 



ISSN: 2249-6645 



i = 2,3,-.. 

7=2,3,.. 



2/ i,„_i,i +4/ u „_ I>2 +4/ n „_ 13 +4/ Un _ 14 + -. + 4/ 11 „_ 1 
/1 w + 2/ ln2 + 2 f { ,„ 3 + 2f, ,„ 4 + ... + 2/, ^ + / IU „ 
2 /i2i 1 + 4/ 1212 + 4/ 12I3 + 4/ 1214 + ... + 4/ 12U _j + 2/ m „ 
4/ 12 2i + 8/ 1222 + 8/ 1223 + 8/ 1224 + ... + 8/ 12 ^ + 4/ m „ 
4/ 12 3i + 8/ 1232 + 8/ I233 + 8/ 1234 + ... + 8/ 123 „_, + 4f l2Jn 



4/,2„-i,, +8/ 12n _ 1>2 +8/ 12 „_ 1>3 +8/ 12n _ 1>4 + ... + 8/ 12n _ ljB _ 1 +4f 12n _ 1>n 
2/,2„i +4/ 12 „ 2 +4/ 12n3 +4/ 12 „ 4 + ... + 4/ 12 „, n _ 1 +2/ 12 „,„ 



/i„n+2/ lnl2 +... + 2/ ln _ 1>B _ 1 +/ lB)n 
2/,„2i + 4/,„ 22 + + 4/ 1 „_ 1>2n _ 1 + 2/ lni . 



2/^_ u +4/ lnn _ li2 + ... + 4/ 1 „,„_ lj „_ 1 +2/ lnB _ 1>n 
/„„ 1 +2/,„, l2 +... + 2/„ + /,„„,„ 

/2ii. + 4/ 2112 + 4/ 2113 + ... + 4/ 2Ull _ 1 +2/ 2U „ 
4/ 2121 + 8/ 2122 + 8/ 2123 + ... + 8/ 212 „_, + 4f 2IZt 



4/ain-i.i +8/ 2 i„-i,2 +8/ 21 „_,, 3 +-- + 8/ 21 „_ 1 ,„_ 1 +4/ 21 „_ ln 
2/ 21 „! +4/ 21 „ 2 +4/ 21 „ 3 + ... + 4/ 21 „„_ 1 +2/ 21 „„ 
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Vim, + 8/2212 + 8/ 2213 + ... + Sf 22ln _, + 4f 22ln 
8/2221 + 16/2222 + 16/ 2223 + ... + 16/ 222n _, + 8f 2 - 

8/223 1 + 1 6/2232 + 16/2233 + - + 16/ 2 23,,-l + 8/ 22 



8/ 22 „_ M +16/ 22M _i, 2 +16/22,^,3 +- + 16/ 2 2„-,,„-, +8/2 
4/ 22b1 +8/ 22 „2 +8/ 22n3 + ... + 8/ 22Bn _, +4/ 22nn 



4/ n -i,„u + 8/„_ 1>nl2 + Sf n _ hnl3 + ... + 8/„_, ,„ .„_! + 4 /„_,„,„ 
8/„_i,„2i + 16/„_ 1 ,„ 22 + 16/^ M + ... + I6f n _ ln2n _, + 8/„_, ,„ , 2 ,„ 



2/„, n - U , + 4/„,„_,, 12 + 4/„ ,„_, >13 + ... + Af nn ,,„ , + 2/,,^ , „ 
4/„,„_ lt21 + 8/ n ,_ lj22 + 8/„.„_ 1 , 23 + ... + S/.,„ ,., , + 4/ Bj „_ 1j2i „ 



4/,..,. , , +8/„,„_ 1 .„_ 1 , 2 +8/„,„_ 1 ,„_ 13 + ... + 8/ n ,„_ 1 ,„_ 1 ,„_ 1 +4/„,„_ 1 ,„_,,„ 
2/,,„-,,„, 1 +4/ Kjn _ ljn , 2 +4/„ jn _ 1> „3 +- + 4/ nK _ lnK _j + 2/ n „_ ln „ 

8 /„-i,„,„-i,i + 16 /„-i,„,„- 1 ,2 +16/„_ 1 .„,„_ 1 ,3 +- + 16/„_ 1 , B ,„_ 1 .„_ 1 +8/„_ ljB> „_ 1>n 
4/„_,, n , n ,i+8/„_ 1;Bi „ >2 +8/„_ 1 „ )H>3 + - + Sf n _ lnnn _ 1 +Sf n _ lnnn 

Ami + 2/„„ 12 + 2/„ nl3 + ... + 2/,„,„_ 1 + /„„,„ 
2/„„2i + 4/„„ 22 + 4/„„23 + ■■■ + 4/ nn2in _j + 4/ nn2 „ 



2/„_ u +4/„_,, 2 +4/„„, M 3 +... + 4/„„ ,„_,,„_, +4f nnn _ l n 

fnnra + 2/„,„ 2 + 2/„„, l3 + ... + 2/„„_ 1 + /„„„,„ 

From these relations the formula (20) is formed, which fully demonstrates the proof of the theorem. 

Now, we calculate the values of integrals by repeated application of quadrature formulas, based on the formula of 

The rectangle formula to calculate the values of the integrals is quite simple. In fact, here<^ ; is taken as 
, X t ] midpoints. For a uniform grid (h i = h, i = 1, n) we obtained formulas of the form 



J 1 =]f(x)d X ~ b ^±f i _ U2 , 

a n 1=1 



(29) 



where f i _ xn = f(x i - —), X Q = a, X n = b . 
Double integral in this case is calculated by the formula 




(30) 



h x - Step on the OX axis, h y y - a step on the OY axis. 

In the same way we define cubature formulas for calculating the values of the three-and four-fold integrals 
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/a = Ijjfix^dxdydz - A ~ aB - b C ~J ±±±f i -u 2 , J -u^-u, > (31) 

^ -f Ijjf^y^dxdydzdu - A ~ a B ~ b ; (32) 

Here /; ^ _ _ ^ = A A,^ _A ^ _ ^ _0. 

■/,-1/2,,-1/a,* 1/2.Z-1/2 J^, 2 2 * 2 , 2 j 



On the basis of approximate formulas for calculating the values of one-, two-, three-and four-fold integrals, 
respectively, expressed by (21) - (24), the formula for calculating the values of n-fold integrals is given . 
Theorem 5. Let the y = f(x l ,X 2 ,...,X n ) function is defined and continuous in a n-dimensional bounded Q„ integration 
domain. Then the cubature formula, obtained by repeated application of the rectangles' formula, has the form 

nj.j.Av.-v. v„./v./v,..,/v„ =n^^zl:...i:/ il _ 1/2 ,,_ 1/2 ,..,„_ 1/2 - 

where x ; e [a t , A i ] , i = 1, ft , 



Proof. The proof is given by induction. 

From (21) - (24) we see that (25) holds for ft = 1;2;3;4 . We assume that (25) is valid for ft = k : 

\\.. D \f{x l ,x 2 ,..., Xk )dx x dx 2 ..Ax k ^ 

Now we show fairness at ft = k + 1 : 

J" J" "' j"/( jc i'- ,c 2 >—>-** ,x k+l )dx l dx 2 ...dx k dx k+l ~ 

The obtained result proofs the theory completely. 

III. CONCLUSIONS 

Thus we have constructed the multi-dimensional cubature formulas to approximate the value of multiple integrals 
by repeated application of the quadrature formula of trapezoids and rectangles to determine the values of multiple integrals, 
which are easy to implement on a computer. 
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